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We report on the calculation of the next-to-leading order QCD corrections to the production of 
top-anti-top quark pairs in association with a hard jet at the Tevatron and at the LHC. We present 
results for the tt+jet cross section and the forward-backward charge asymmetry. The corrections 
stabilize the leading-order prediction for the cross section. The charge asymmetry receives large 
corrections. 
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INTRODUCTION 

More than ten years after its discovery, the dynamics 
and many properties of the top quark, such as its elec- 
troweak quantum numbers, are not yet precisely mea- 
sured. Since the top quark is the heaviest of the known 
elementary particles, it is widely believed that it plays a 
key role in extensions of the Standard Model. This ren- 
ders experimental investigations of the top quark partic- 
ularly important. Up to now the main (direct) source of 
information on top quarks are top-quark pairs produced 
at the Tevatron. Only recently first evidence for single- 
top production has been found It is important to 
note that in the inclusive tt sample a significant fraction 
comprises tt+jet events. An investigation of the process 
of tt production in association with a hard jet can, thus, 
certainly improve our knowledge about the top quark. 

In this context, the forward-backward charge asym- 
metry of the top (or anti-top) quark 0, H, 0, @ is of 
particular interest. In inclusive tt production it appears 
first at one loop, because it results from interferences of 
C-odd with C-even parts of double-gluon exchange be- 
tween initial and final states. This means that the avail- 
able prediction for tt production — although of one-loop 
order — describes this asymmetry only at leading-order 
(LO) accuracy. In tt+jet production the asymmetry 
appears already in LO. Thus, the next-to-leading order 
(NLO) calculation described in the following provides a 
true NLO prediction for the asymmetry. Our calculation 
will, therefore, be an important tool in the experimen- 
tal analysis of this observable at the Tevatron where the 
asymmetry is measureable as discussed in Ref. @. 

Measuring the cross section of the related process of 
tt+7 production provides direct access to the electric 
charge of the top quark. Obviously NLO QCD predic- 
tions to this process are important for a reliable analysis. 
They can be obtained from tt+jet production presented 



in this letter via simple substitutions. We postpone this, 
however, to a forthcoming publication. 

Finally, a signature of tt in association with a hard jet 
represents an important background process for searches 
at the LHC, such as the search for the Higgs boson in the 
weak-vector-boson fusion or ttH channels. 

The above-mentioned issues clearly underline the case 
for an NLO calculation for tt+jet production at hadron 
colliders. We report here on a first calculation of this 
kind; a status report was already given in Ref. 0. 



DETAILS OF THE NLO CALCULATION 

At LO, hadronic tt+jet production receives contribu- 
tions from the partonic processes qq — > tig, qg — > ttq, 
qg — > tiq and gg — > tig. The first three channels are re- 
lated by crossing symmetry to the amplitude — > tiqqg. 

Evaluating 2^3 particle processes at the NLO level, 
is non-trivial, both in the analytical and numerical parts 
of the calculation. In order to prove the correctness of 
our results we have evaluated each ingredient twice using 
independent calculations based — as far as possible — on 
different methods, yielding results in mutual agreement. 



Virtual corrections 

The virtual corrections modify the partonic processes 
that are already present at LO. At NLO these correc- 
tions are induced by self-energy, vertex, box (4-point), 
and pentagon (5-point) corrections. The prototypes of 
the pentagon graphs, which are the most complicated 
diagrams, are shown in Figure [1] 

Version 1 of the virtual corrections is essentially ob- 
tained following the method described in Ref. Q, where 
ttH production at hadron colliders was considered. Feyn- 
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FIG. 1: A representative set of pentagon diagrams for 
hadronic tt+jet production at NLO QCD. 



man diagrams and amplitudes have been generated with 
the FeynArts package and further processed with in- 
house Mathematica routines, which automatically create 
an output in Fortran. The IR (soft and collinear) singu- 
larities are analytically separated from the finite remain- 
der as described in Refs. @, E3]- The tensor integrals 
appearing in the pentagon diagrams are directly reduced 
to box integrals following Ref. 



This method does 



not introduce inverse Gram determinants in this step, 
thereby avoiding notorious numerical instabilities in re- 
gions where these determinants become small. Box and 
lower-point integrals are reduced a la Passarino-Veltman 
[l2T ] to scalar integrals, which are either calculated ana- 



lytically or using the results of Refs. [HI, 03, EH. Suffi- 
cient numerical stability is already achieved in this way. 
Nevertheless the integral evaluation is currently further 
refined by employing the more sophisticated methods de- 
scribed in Ref. 16] in order to numerically stabilize the 
tensor integrals in exceptional phase-space regions. 

Version 2 of the evaluation of loop diagrams starts 
with the generation of diagrams and amplitudes via 
QGRAF 17], which are then further manipulated with 
Form (lij and eventually automatically translated into 
C++ code. The reduction of the the 5-point tensor in- 
tegrals to scalar integrals is performed with an extension 
of the method described in Ref. [l9j| ■ In this procedure 
also inverse Gram determinents of four four-momenta are 
avoided. The lower-point tensor integrals are reduced 
using an independent implementation of the Passarino- 
Veltman procedure. The IR-finite scalar integrals are 
evaluated using the FF package 2^, 21 1. Although the 
entire procedure is sufficiently stable, further numerical 
stabilization of the tensor reduction is planned follow- 
ing the expansion techniques suggested in Ref. [22| for 
exceptional phase-space regions. 



Real corrections 

The matrix elements for the real corrections are given 
by — > itgggg, — > ttqqgg, — > ttqqq'q' and — > ttqqqq. 
The various partonic processes are obtained from these 
matrix elements by all possible crossings of light particles 
into the initial state. 

The evaluation of the real-emission amplitudes is per- 
formed in two independent ways. Both evaluations em- 
ploy (independent implementations of) the dipole sub- 
traction formalism 
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|25| for the extraction of IR 
singularities and for their combination with the virtual 
corrections. 

Version 1 results from a fully automated calculation 
based on helicity amplitudes, as described in Ref. [261 ]. 
Individual helicity amplitudes are computed with the 
help of Berends-Giele recurrence relations j2?| ■ The eval- 
uation of color factors and the generation of subtraction 
terms is automated. For the channel gg — ► tigg a dedi- 
cated soft-insertion routine [28| is used for the generation 
of the phase space. 

Version 2 uses for the LO 2^3 processes and the 
gg — ► tigg process optimized code obtained from a Feyn- 
man diagramatic approach. As in version 1 standard 
techniques like color decomposition and the use of helic- 
ity amplitudes are employed. For the 2 — > 4 processes in- 
cluding light quarks, Madgraph [2£| has been used. The 



subtraction terms according to Ref. [25| are obtained in 
a semi-automatized manner based on an in-house library 
written in C++. 



NUMERICAL RESULTS 



In the following we consistently use the CTEQ6 13, H3 
set of parton distribution functions (PDFs). In detail, 
we take CTEQ6L1 PDFs with a 1-loop running a s in 
LO and CTEQ6M PDFs with a 2-loop running a s in 
NLO. The number of active flavours is Np — 5, and 
the respective QCD parameters are Aj ° = 165 MeV and 
A^ s = 226 MeV. Note that the top-quark loop in the 
gluon self-energy is subtracted at zero momentum. In 
this scheme the running of a s is generated solely by the 
contributions of the light quark and gluon loops. The 
top-quark mass is renormalized in the on-shell scheme, 
as numerical value we take m t = 174 GeV. 

We apply the jet algorithm of Ref. [32| with R = 1 for 
the definition of the tagged hard jet and require a trans- 
verse momentum of PTjct > 20 GeV for the hardest jet. 
The outgoing (anti-)top quarks are neither affected by 
the jet algorithm nor by the phase-space cut. Note that 
the LO prediction and the virtual corrections are not in- 
fluenced by the jet algorithm, but the real corrections are. 

Figure [2] shows the dependence of the integrated LO 
cross sections on the renormalization and factorization 
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FIG. 2: Scale dependence of the LO cross sections for tt+jet 
production at the Tevatron (upper plot) and at the LHC 
(lower plot), where the renormalization and factorization 
scales are set equal to fi. The individual contributions of 
the various partonic channels are shown also separately. 
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FIG. 3: Scale dependence of the LO and NLO cross sections 
for tt+jet production at the Tevatron (upper plot) and at the 
LHC (lower plot), where the renormalization and factoriza- 
tion scales are set equal to ji. 



scales, which are identified here (/i = /i rcn = /if ac t)- The 
variation ranges from /i = 0.1 m t to (i = 10 m t . The 
dependence is rather large, illustrating the well-known 
fact that the LO predictions can only provide a rough 
estimate. At the Tevatron the qq channel dominates the 
total pp cross section by about 85%, followed by the gg 
channel with about 7%. At the LHC, the gg channel 
comprises about 70%, followed by qg with about 22%. 

In Figure [3] the scale dependence of the NLO cross 
sections is shown. For comparison, the LO results are 
included as well. As expected, the NLO corrections sig- 
nificantly reduce the scale dependence. We observe that 
arround (i w m t the NLO corrections are of moderate 
size for the chosen setup. 

Finally, we have performed a first study of the forward- 
backward charge asymmetry of the top quark at the 
Tevatron. In LO the asymmetry is defined by 
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FB.LO 
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= cr LO (y t >0)±erLo(2/t<0), (1) 



where y t denotes the rapidity of the top quark. Cross- 
section contributions a(yt < 0) correspond to top quarks 
in the forward or backward hemispheres, respectively, 
where incoming protons fly into the forward direction 
by definition. Denoting the corresponding NLO contri- 
butions to the cross sections by Sa^ LO , wc define the 
asymmetry at NLO by 



^FB.NLO 



'LO 



'LO 



1 



5a 



NLO 



'LO 



^NLO 



'LO 



(2) 



i.e. via a consistent expansion in a s . Note, however, that 
the LO cross sections in Eq. ([2]) are evaluated in the NLO 
setup (PDFs, a s ). 

Figure U shows the scale dependence of the asymmetry 
at LO and NLO. The LO result for the asymmetry is of 
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FIG. 4: Scale dependence of the LO and NLO forward- 
backward charge asymmetry of the top quark in pp — > 
tt+jet+X at the Tevatron, where the renormalization and 
factorization scales are set equal to fi. 

order oP s and does therefore not depend on the renormal- 
ization scale. The plot for the LO result shows a mild 
residual dependence on the factorization scale, but the 
size of this variation does not reflect the theoretical un- 
certainty, which is much larger. The NLO corrections 
to the asymmetry are of order a l s and depend on the 
renormalization scale. It is therefore natural to expect 
a stronger scale dependence of the asymmetry at NLO 
than at LO, as seen in the plot. The size of the asym- 
metry, which is about —7% at LO, is drastically reduced 
by the NLO corrections. The residual scale dependence 
suggests that the asymmetry is — 1.5(±1.5)% for the cho- 
sen setup. Before closing, we comment on an alternative 
possibility to define the asymmetry via the unexpanded 
ratio of the NLO-corrected cross sections Clo + ^ a NLO ■ 
With this definition (not shown in the plot), one observes 
an even larger scale dependence with a dramatic rise for 
low scales fi where the NLO cross section (the denomina- 
tor in the asymmetry) becomes perturbatively unstable 
and turns to negative values (see Figured]). 

CONCLUSIONS 

First predictions for tt+jet production at hadron col- 
liders have been presented at NLO QCD. For the cross 
section the NLO corrections drastically reduce the scale 
dependence of the LO predictions, which is of the order 
of 100%. The charge asymmtry of the top quarks, which 
is going to be measured at the Tevatron, is significantly 
decreased at NLO and is almost washed out by the resid- 
ual scale dependence. 
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